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ABSTRACT 


This  report  describes  some  of  the  algorithms  used  in  a  prototype  automatic  binary 
image  scaling  facility.  This  facility  was  designed  and  implemented  to  encode  scanned, 
digitized,  binary  images  and  to  generate  faithful  reproductions  of  the  images  in  different 
sizes  and  at  different  resolutions. 


CHAPTER  1.  INTRODUCTION 

This  report  describes  some  of  the  basic  algorithms  for  automatic  scaling  of  scanned 
binary  images.  Typically,  a  binary  image  may  represent  a  digitized  font  or  logo.  It  is 
scanned  into  a  text  formatting  or  typesetting  system  and  is  later  displayed  or  printed, 
frequently  in  different  point  sizes  and  on  devices  having  different  resolutions.  The 
algorithms  described  here  are  used  in  a  prototype  automatic  scaling  facility  (ASF) 
designed  and  implemented  to  carry  out  the  following  tasks: 

(1)  Given  the  bit-map  representation  of  a  binary  image,  detects  its  boundary  and 
generates  an  efficient  encoding  of  the  boundary. 

(2)  Finds  a  set  of  piecewise  polynomials  which  approximates  the  boundary  with 
specified  error  tolerance. 

(3)  Scales  the  approximate  boundary  to  a  reduced  size. 

(4)  Generates  a  bit-map  representation  of  the  image  at  the  desired  resolution  by 
inking  the  interior  of  the  boundary  produced  in  the  step  above.  This 
transformed  bit-map  representation  may  be  used  as  input  to  a  printer  or 
display. 

The  prototype  ASF  was  designed  to  be  used  with  a  limited  amount  of  user  intervention. 
Specifically,  the  tasks  of  choosing  the  first  point  on  the  image  boundary,  marking  up 
unconnected  portions  of  the  same  image,  and  judging  the  visual  quality  of  an 
approximate  boundary  are  carried  out  by  the  user.  Appendix  A  describes  an  example 
illustrating  a  way  in  which  this  facility  may  be  used. 

Usually,  fewer  bits  are  required  to  encode  and  store  a  binary  image  when  it  is 
specified  by  its  boundary  points  or  by  piecewise  polynomials  approximating  the 
boundary.  Hence,  the  technique  described  here  can  also  be  used  to  achieve  data 
compression  in  representations  of  binary  images.  However,  our  primary  objective  is  to 
generate  faithful  reproductions  of  the  binary  images  in  different  sizes  and  at  different 
resolutions.  Minimization  of  the  storage  requirement  is  not  an  important  consideration 
here. 


The  problem  of  scaling  binary  images  to  different  sizes  or  translating  them  to 
different  resolutions  has  been  studied  using  various  techniques  in  digital  filtering  and  in 
contour  smoothing  and  adjustment  [1,  2].  In  these  studies,  image  resolution  is  translated 
by  either  insertion  or  deletion  of  pixels  and  by  smoothing  the  contours  of  the  resultant 
image.  The  algorithms  described  in  this  report  were  designed  to  achieve  the  same 
objectives  as  in  these  two  studies.  The  approach  taken  here  is  to  first  find  the  boundary 
of  the  given  scanned  image,  represent  the  boundary  by  piecewise  polynomials,  reduce  the 
boundary  by  scaling  the  polynomials,  and  then  ink  the  interior  of  the  boundary  at  the 
desired  resolution. 

The  use  of  piecewise  polynomials  to  represent  the  boundary  of  a  binary  image  was 
suggested  in  TEX  [3].  In  the  case  of  TEX,  mathematical  representation  of  font 
boundaries  provides  a  font  designer  the  needed  tool  for  shaping  and  adjusting  the  form 
of  the  character  until  a  desirable  form  is  found.  In  our  case,  the  piecewise  polynomials 
are  used  to  approximate  given  boundaries.  Because  of  this  difference  in  our  objectives, 
different  polynomials  were  found  to  be  more  suitable  for  boundary  representation.  The 
relative  merits  of  using  different  piecewise  polynomials  to  approximate  a  curved 
boundary  have  been  evaluated  extensively.  References  [4  -  7]  provide  the  bases  for  our 
choices  of  polynomials. 

In  Chapter  2,  the  algorithm  for  finding  the  boundary  of  a  binary  image  is 
described.  The  types  of  polynomials  used  to  approximate  image  boundaries  are  described 
in  Chapter  3  together  with  the  justifications  for  choosing  them.  The  error  measure  used 
to  rank  the  merits  of  approximate  boundaries  and  the  equations  relating  the  parameters 
of  the  polynomials  to  coordinates  of  boundary  points  are  also  described  in  this  chapter. 
The  descriptions  of  the  equipment  on  which  the  prototype  ASF  was  implemented  and 
the  programs  used  in  this  facility  can  be  found  in  Appendices  A  and  B.  Chapter  4 
contains  a  summary  of  the  work  described  in  this  report.  Several  illustrative  examples 
are  given  in  this  chapter. 
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CHAPTER  2.   BOUNDARY  DEFINITION 

This  chapter  describes  the  algorithm  used  to  find  the  boundary  of  a  binary  image 
from  the  bit-map  representation  of  the  image.  In  particular,  the  algorithm  is  applied  to 
each  separate  component  of  the  binary  image.  When  the  given  binary  image  has  more 
than  one  component,  (e.g.,  the  letter  "•"  "),  the  boundary  of  each  component  is  found 
independently.  Similarly,  in  the  case  where  the  white  region  is  not  connected,  (e.g,  the 
letter  "O"),  the  two  disjoint  boundaries  are  found  separately.  The  relative  positions  of 
the  disjoint  boundaries  and  individual  components  are  specified  by  the  coordinates  of 
one  point  on  each  of  the  boundaries.  In  this  chapter,  our  attention  is  confined  to  the 
process  of  finding  one  boundary  when  a  point  on  the  boundary  is  given. 

We  view  each  pixel  as  a  small  square.  The  coordinates  of  the  pixel  are  the 
coordinates  of  the  center  of  the  square.  For  a  given  bit-map  representation  of  a  binary 
image,  the  value  of  a  square  is  1  if  it  is  part  of  the  binary  image.  The  value  of  a  square 
within  the  white  space  around  the  image  is  0.  Let  #,•  be  a  square  whose  value  is  equal  to 
1  and  which  shares  at  least  one  side  with  another  square  whose  value  is  equal  to  0.  An 
arrow  drawn  along  this  common  side  is  called  a  boundary  vector  associated  with  the  5,. 
(For  the  sake  of  concreteness,  the  direction  of  the  boundary  vector  is  always  chosen  so 
that  as  one  transverses  along  and  in  the  direction  of  the  boundary  vector,  the  square 
with  value  1  is  on  the  left  of  the  path.)  In  the  example  shown  in  Figure  1,  each  of  the 
squares  Blt  B2,  and  Z?s  shares  at  least  one  side  with  a  square  whose  value  is  0.  The 
boundary  vectors  associated  with  these  squares  are  as  shown. 

Two  squares  are  said  to  be  adjacent  to  one  another  when  they  shares  one  side  or  a 

corner.    Let  Bit  for  t  =  1,  2,  ,  N,  be  N  squares  each  of  which  shares  one  side  with  a 

square  whose  value  is  0.  Moreover,  these  N  squares  are  adjacent  to  each  other.  We  call 
these  squares  the  boundary  points  of  the  binary  image  if,  as  one  transverses  along  a  path 

formed  by  the  boundary  vectors  associated  with  the  squares  Bx,  B2,  ,  and  BN,  the 

squares  on  the  left  of  the  path  have  value  1  and  the  squares  on  the  right  of  the  path 
have  value  0.  For  the  images  shown  in  Figure  1,  squares  with  value  1  are  shaded. 
Boundary  points  are  darker  shaded  squares.  The  closed  path  formed  by  the  connecting 
the  centers  of  the  boundary  points  is  called  the  boundary  of  the  binary  image. 

Given  a  binary  image,  our  problem  is  to  find  its  boundary.  The  boundary  finding 
algorithm    used   in   the   prototype   ASF,   referred   to   as   Algorithm   BF,   works   in   the 
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Figure  1.   Examples  of  Boundary  Points  and  Boundary  Vectors 


following  manner: 

Algorithm  BF.   The  boundary  points  found  by  the  algorithm  are  kept  in  the  boundary 
stack.  The  algorithm  works  as  follows: 

•  The  algorithm  begins  with  a  square  that  is  known  to  be  a  boundary  point.  Also 
given  is  a  boundary  vector  associated  with  the  boundary  point.  They  are 
referred  to  as  the  current  boundary  point  and  current  boundary  vector, 
respectively. 

•  In  order  to  find  the  next  boundary  point,  four  situations  shown  in  Figure  2  are 
considered,  (a).  In  Case  a,  the  three  squares  adjacent  to  the  current  boundary 
point  have  the  value  0  and  the  current  boundary  vector  is  as  shown  on  the 
left.  The  next  step  is  to  change  the  current  boundary  vector  to  be  the  one 
shown  on  the  right.  No  new  point  is  added  to  the  boundary  point  stack,  (b). 
In  Case  b,  the  current  boundary  vector  is  changed  as  shown.  Square  No.  1  is 
added  to  the  boundary  point  stack,  (c)  and  (d).  In  the  Cases  of  c  and  d,  the 
current  boundary  vector  is  changed  as  shown.  Square  No.  2  is  added  to  the 
boundary  point  stack. 

•  The  algorithm  terminates  when  the  starting  boundary  vector  becomes  the 
current  boundary  vector  again. 

Algorithm  BF  is  conceptually  simple.  It  is  implemented  by  the  subroutine  called  the 
Boundary  Eater  listed  in  Appendix  B.  This  program  is  complex  (see  lines  10000 
through  10700)  due  to  the  fact  that  all  situations  symmetrical  to  the  four  cases  shown  in 
Figure  2  must  be  dealt  with.  For  example:  the  situation  shown  in  Figure  3  is 
symmetric  to  Case  d  illustrated  in  Figure  2.  To  deal  with  all  of  the  cases,  a  concept  of 
a  greater  or  a  smaller  relationship  between  a  boundary  point  and  its  associated  boundary 
vector  is  introduced.  When  the  boundary  point  is  below  or  is  to  the  right  of  the  current 
boundary  vector,  the  point  is  said  to  be  greater  than  the  boundary  vector.  When  the 
boundary  point  is  above  or  is  to  the  left  of  the  boundary  vector,  the  point  is  said  to  be 
smaller  than  the  vector.  Furthermore,  with  the  current  boundary  point  being  square  0  in 
Figure  2,  square  2  is  referred  to  as  the  diagonally  opposite  point.  Square  1  is  referred  to 
as  the  other  point  of  interest.  Thus,  the  exact  action  taken  by  Boundary  Eater  depends 
on: 
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Figure  2.   Four  Situations 
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Figure  3.   A  Case  symmetrical  to  Case  d. 


(1)  The  direction  of  the  boundary  vector  (stored  in  variable  DR  as  U,  D,  L,  or  R 
for  up,  down,  left,  or  right,  respectively,) 

(2)  The  value  of  the  diagonally  opposite  point  (stored  in  variable  DO  as  1  or  0,) 

(3)  The  value  of  the  other  point  of  interest  (stored  in  variable  OT  as  1  or  0,) 

(4)  The  relationship  between  the  current  boundary  point  and  the  boundary  vector 
(stored  in  variable  RL  as  either  G  or  S  for  greater  or  smaller,  respectively.) 

Observe  that  the  coordinates  of  the  diagonally  opposite  point  and  the  other  point  of 
interest  are  functions  of  the  coordinates  of  the  current  boundary  point,  the  direction 
DR    of  the  current  boundary  vector,  and  the  relationship  RL . 

Specifically,  the  subroutine  Boundary  Eater  requires  the  following  data  as  inputs: 
(1)  the  coordinates  of  the  starting  boundary  point,  (X(0),  Y(0)),  (2)  the  direction  of  the 
starting  boundary  vector,  U,D ,L,  or  R ,  and  (3)  the  relationship  between  the  starting 
boundary  point  and  boundary  vector,  G  or  S.  The  subroutine  returns  the  coordinates 
of  the  boundary  points  in  arrays  X  and  Y  from  the  starting  point  given  by  ( X(0),  Y(Q) ) 
to  the  ending  point  given  by  (X(TP),  Y(TP))  (which  is  the  same  as  the  starting  point.) 
Whether  the  boundary  defines  a  clockwise  or  counter-clockwise  contour  is  reflected  by 
the  coordinates  of  adjacent  points,  (X(K  +  1),  Y(K  +  1))  and  ( X{K),  Y(K)).  The 
coordinates  of  the  boundary  points  in  arrays  X  and  Y  are  subsequently  used  as  inputs 
to  curve  fitting  algorithms  discussed  in  the  next  chapter. 

Since  adjacent  array  elements  contain  coordinates  of  adjacent  boundary  points,  it  is 
only  necessary  to  store  the  full  coordinates  of  the  first  boundary  point.  Given  the 
coordinates  of  a  boundary  point,  which  one  of  the  8  squares  adjacent  to  this  point  is  the 
next  boundary  point  may  be  specified  using  3  bits.  This  technique  for  encoding 
boundaries  of  binary  images  is  known  as  chained  link  encoding  and  has  been  used  to 
achieve  data  compression.  For  example,  the  image  in  Figure  4  is  roughly  130  squares 
high  by  200  squares  wide,  requiring  26,000  bits  if  stored  as  a  bit  map.  However,  the 
boundary  consists  of  only  540  squares.  If  the  boundary  is  stored  using  3  bits/square,  only 
1,620  bits  are  required  to  store  this  image. 


Figure  4.   An  Example 


CHAPTER  3.  ALGORITHMS  TO  FIND  APPROXIMATE  BOUNDARY 

This  chapter  describes  the  algorithm  used  to  find  an  approximate  boundary  from 
the  given  coordinates  of  boundary  points.  In  particular,  the  choices  of  the  piecewise 
polynomials  used  to  specify  the  approximate  boundary  and  the  error  measure  used  to 
judge  the  quality  of  the  approximation  are  justified.  The  algorithm,  referred  to  as  the 
Boundary  Compacter,  was  designed  to  be  used  in  an  interactive  environment.  After  the 
boundary  of  a  binary  image  is  found,  the  user  would  specify  the  amount  of  maximum 
tolerated  error  in  the  approximate  boundary  and  call  the  Boundary  Compacter  to  find 
an  approximate  boundary  with  error  less  than  the  specified  maximum  tolerated  error. 
After  comparing  the  resultant  approximate  boundary  with  the  original,  he  may  wish  to 
repeat  the  approximation  process  with  a  smaller  maximum  tolerated  error. 

Figure  5  illustrates  a  frequently  encountered  example  of  boundary  points  defining  a 
segment  of  the  image  boundary.  Note  that  the  boundary  points  lie  on  a  curve  which  is 
not  a  function  of  X  or  of  Y  .  One  approach  to  finding  the  piecewise  polynomial 
approximation  of  such  a  curve  is  to  segment  the  curve  so  that  some  pieces  are  functions 
of  X  and  some  pieces  are  functions  of  Y .  However,  this  approach  may  force  one  to 
choose  unnecessarily  short  segments.  Moreover,  a  great  deal  of  testing  is  required  in 
order  to  determine  which  segments  are  functions  of  X  and  which  are  functions  of  Y. 

The  approach  taken  here  is  as  follows:  Since  different  pieces  of  polynomials  are  used 
to  approximate  segments  of  the  boundary,  the  straight  line  connecting  the  endpoints  of 
each  segment  may  be  used  to  establish  a  new  X-axis.  That  segment  of  the  boundary 
may  be  approximated  by  a  polynomial  of  this  new  variable  X.  This  approach  is 
illustrated  in  Figure  6.  We  note  that  for  each  pair  of  endpoints,  (X$  Y$)  and  (Xg  Yg), 
defining  a  segment  of  the  boundary,  there  corresponds  a  coordinate  system  defined  as 
follows:  The  new  X  axis  is  the  line  joining  the  two  endpoints  and  the  new  point  of  origin 
is  the  point  (Xs  Ys).  The  coordinates  of  each  point,  (XK  YK),  on  the  boundary  become 
the  coordinates  (XB tYB).  We  can  now  proceed  to  find  the  coefficients  of  the  polynomial 
P(XB)  so  that 

(1).  P(0)  =  0  and  P(XB)  =  0,  and 
(2).  Z[P{XB)-  YB  ]2  is  minimized. 
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Figure  5.  An  Example  of  Boundary  Segment 


(XsY. 


Figure  6.  Definition  of  the  New  Coordinate  System 
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3-1.  Choice  of  the  Polynomial  P(XB) 

The  polynomials  chosen  to  approximate  segments  of  the  image  boundary  are 
parabolic  segments  which  can  be  oriented  as  desired.  Such  a  parabola  is  given 
parametrically  by  the  following  two  equations. 

X{t)  =  at{t-1)+  0t  +  Xs 

Y(t)  =  1t(t-1)+  6t+  Ys 

where  /?=  Xg  -Xs  and  6=  YE  -  Y$.  The  coefficients  a  and  7  can  be  determined  with  a 
little  more  work:  Let 

Ts=Vsp-Us6 

TE  =  VE0-UE6 

Dj^VsUE-VEUs 
Then 

a  =  (t/£rs+  USTE)/Dj 

1  =  (VETS+  VSTE)I  Dj 

The  last  two  equations  will  not  yield  finite  solutions  when  Dj  =  0  which  occurs  only 
when  VSUE  =  VE  Us  .  That  is,  when  the  tangent  lines  at  the  two  endpoints  of  the 
boundary  segment  are  parallel.  Note  that  the  subroutine  Quardratic  in  Appendix  E 
does  not  deal  with  this  situation.  If  Dj  =  0  ,  it  merely  stops.  (See  line  8080  of  the 
program.) 

A  variety  of  smooth  curves  may  be  created  by  chaining  segments  of  parabolas 
together  with  matching  slopes  at  common  endpoints.  On  the  other  hand,  a  sharp  corner 
at  some  point  can  be  created  by  joining  two  parabolas  at  a  common  endpoint  with 
different  slopes  at  that  point.  Indeed,  it  is  more  advantageous  for  us  to  use  parabolas 
here  instead  of  cubic  polynomials  of  the  types  suggested  in  [3  -  6].  Such  a  cubic 
polynomial  may  contain  an  inflection  point  which  may  be  difficult  to  eliminate.  It  is 
relatively  easy  to  approximate  a  segment  with  an  inflection  point  by  using  two  parabolas 
as  shown  in  Figure  6.  Since  minimizing  storage  is  not  one  of  our  objectives,  the  need  for 
one  more  polynomial  in  this  case  is  not  an  objection. 
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Alternatively,  we  may  use  cubic  polynomials  with  the  general  form 

P(X)  =  (AX+  B)X(X-Bb)  (1) 

where  Bb  is  as  defined  in  Figure  6.  For  notational  simplicity,  X  is  used  to  mean  Xq 
hereafter  unless  it  is  stated  otherwise.  Similarly,  whenever  the  coordinates  XK  and  YK 
are  referred  to  hereafter,  we  do  not  mean  the  original  coordinates  of  the  boundary  point. 
Rather,  they  refer  to  the  coordinates  with  respect  to  the  new  coordinate  system  defined 
in  Figure  6.  (Lines  18050  and  18060  of  the  program  listed  in  Appendix  B  show  how  to 
compute  the  new  coordinates  from  the  original  coordinates.)  Note  that  in  the  range  of 
(0,  Bb),  P(X)  has  no  inflection  point. 

Hereafter,  we  confine  our  discussions  to  the  cubic  polynomials  of  the  form  given  by 
Eq.(l)  unless  stated  otherwise.  Similar  statements  and  conclusions  can  be  made  for  the 
case  when  the  quardratic  polynomials  are  used. 

3-2.  Determination  of  the  Polynomial  Coefficients 

Given  P(0)  ===  0  and  P{Bb)  =  0,  we  find  the  polynomial  of  the  form  given  by  Eq. 
(1)  which  interpolates  {XSYS)  and  (XE,YE).   Let 

g(A,B)  -  UP(XK)-  Yk]2  =  U(AXK  +  B)XK{XK-Bby  YK]2 

K  K 

We  wish  to  determine  the  values  of  A  and  B  which  minimize  g(A  ,B).  Setting  the 
partial  derivatives 

-l3-=Yl2[{AXK+  B)XK{XK-Bb)-YK}XK\XK-Bb) 
dA         K 

J£-=Yt2\{AXK+B)XK{XK-Bb)-YK\XK{XK-Bb) 
db         K 

to  zero,  we  obtain 

A  E  XK*(XK-Bb)2+  BY,  XK*(XK-Bb)2  =  £  YKXK\XK-Bb) 

A  EV(4-B»)2+  BYlXK2{XK-Bbf  =  £  YKXK{XK-Bb) 

Rewriting  these  two  equations  using  the  notation  in  the  subroutine  Least  Squares  Cubic 
Coeff,  we  have 

A  Sb  +    BSi  =  S2  (2a) 


13 


A  S4  +    B  St  =  St  (2b) 

These  two  equations  can  easily  be  solved  for  A  and  B.  (See  lines  18160  -  18180  of  the 
program  in  Appendix  B.) 

3-3.  Finding  A  Better  Approximate  Boundary 

Once  a  polynomial  segment  is  found  to  approximate  a  boundary  segment,  it  is 
necessary  to  determine  how  good  the  approximation  is.  For  the  sake  of  tractability,  the 
polynomial  is  chosen  here  to  minimize  the  sum  of  the  square  errors  Yj[^(^k)~  Yk]2  • 
However,  this  error  measure  is  not  a  good  criterion  when  used  to  judge  the  "nearness"  of 
B(X)  to  the  actual  boundary.    Rather,  the  maximum  absolute  error,  max  \P(XK)-  YK\ 

should  be  used.  Indeed,  the  maximum  absolute  error  is  the  measure  used  in  the 
subroutine  Average  Absolute  and  Maximum  Error.  The  value  of  the  maximum  absolute 
error  in  the  polynomial  found  to  approximate  a  boundary  segment  is  returned  by  the 
subroutine  in  the  variable  UR  .  (This  subroutine  re-calculates  the  values  of  the  new 
coordinates  from  the  original  coordinates  of  every  boundary  point  because  the  values  of 

the  new  coordinates  were  not  saved  by  the  Least  Squares  Cubic  Coeff  Subroutine 

see  lines  18440  and  18450.   The  value  of  P{XK)  is  given  by  YT  in  the  line  18400.) 

In  order  to  allow  for  later  drawing  of  the  polynomial,  the  coordinates  of  the 
starting  point,  (XS,YS),  the  coordinates  of  the  endpoint,  (XE,  YE),  and  the  coefficients 
A  and  B  are  retained.  The  segment  of  the  polynomial  is  drawn  using  the  following 
program  segment: 

Bl  =  X(E)-X{S):  B2=Y(E)-Y(S) 
For  XT  =  0  to  Bb     Step  H 

YT=(A  *DT+  B)XT  *(XT-Bb) 

Pl  =  X(S)+  (YT  *B2  +  XT  *Bi)/Bb 

P2=Y{S)  +  {XT  *B2-  YT  *Bi)j  Bb 

plot '(jyy 

Next  XT 

The  subroutine  Better  Plotter  works  in  a  similar  manner. 

We  note  that  the  size  of  the  step  H  must  be  determined  to  give  an  unbroken  plot. 
To  choose  a  reasonable  value  of  H,  we  consider  the  general  problem  of  plotting  the  point 
between  two  points  (Xl,Yl)    and  (X2,Y2).    Roughly,  one  is  inclined  to  do  something 
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such  as 

For  T  =  0   to  1      Step// 

PlotX(r),  Y{T) 
Next  T 

The  value  of  H  must  be  selected  so  that  the  resulting  plot  is  continuous  on  an  integer 
grid.  If  the  point  (X(T),  Y(T})  has  already  been  plotted,  then  the  next  point  to  be 
plotted  will  be  X(T  +  //),  Y(T  +  //).  To  ensure  a  "continuous"  plot  on  an  integer  grid, 
we  must  have  \X(t  +  H)-X(t)\  <  1  and  |  Y{t  +  H)-Y(l)\  <  1  .  In  other  words,  if 
the  parabolas  of  the  form  defined  in  Section  3-2  are  used, 

\X{t+H)-X(t)\  =  \X*(e)\  \H\ 

by  the  mean  value  theorem,  where  c  E{t,t  +  H)  <  (0,1),  we  have 

\X(t+  H)-X(t)\  =  \X'(e)\  \H\  =  \2ae+  (/?-  alpha)  |  |//|. 

The  right  hand  side  is  less  than  or  equal  to 

max[|/?-a|  \H\,\0+a\  \H\]. 

for  eG  (0,1).    Similarly, 

\Y(t+H)-Y(t)\<mBx[\6-1\  |H|,|*+7|  \H\]. 

Therefore,  to  ensure  that  |X(f  +  H)-X(t)\  <  1  and  |  Y{t  +  H)-  Y{t)\  <  1,  we  should 
select  H  so  that  H  =  1  /  MX  where 

MY  =  max[|0-a|,|/?+a|,|*-7|,|*+'y|]. 

However,  this  choice  of  H  is  overly  pessimistic  resulting  in  too  many  points  plotted.  For 
this  reason,  in  lines  8130  -  8260  of  the  subroutine  Better  Plotter,  where  the  value  of  MX 
is  calculated  for  cubic  polynomials  given  by  Eq.  (1).  H  is  chosen  to  be  //  =  1.9/  MX  . 
This  value  was  arrived  at  experimentally.  To  guard  against  occasional  gaps,  the  rest  of 
the  plotting  subroutine  does  not  plot  individual  points  but  connects  successive  points 
with  short  line  segments. 

3-4.  Boundary  Compacter 

Boundary  Compacter  is  the  program  which  calls  the  other  subroutines  to  create  a 
sequence  of  polynomial  segments.  Together,  these  segments  approximate  the  entire 
boundary  with  maximum  absolute  error  less  than  the  user  specified  maximum  tolerated 
error. 
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The  input  to  the  Boundary  Compacter  is  the  boundary  point  stack  X  and  Y 
generated  by  the  Boundary  Eater.  The  Boundary  Compacter  starts  with  the  first  four 
boundary  points.  It  determines  the  parameters  A  and  B  according  to  Eq.  (2).  These 
parameters  specify  the  cubic  polynomial  segment  used  to  approximate  the  boundary 
defined  by  these  four  points.  Then  the  maximum  absolute  error,  UR  ,  is  calculated.  If 
UR  is  less  than  the  given  maximum  tolerated  error,  denoted  by  UP  ,  one  more  point  is 
adjoined  to  the  original  four  and  the  entire  process  is  repeated.  In  this  way,  the 
polynomial  segment  grows  until  the  inclusion  of  an  additional  boundary  point  causes 
UR  to  exceed  UP.  The  starting  point  and  the  endpoint  of  the  cubic  segment  and  the 
parameters.  A  and  B,  defining  the  cubic  polynomial  are  saved  so  that  the  segment  can 
be  drawn  later.  Then  a  new  cubic  segment  is  grown  starting  with  four  boundary  points 
not  previously  considered. 

Growing  an  approximate  boundary  segment  a  point  at  a  time  may  be  very  slow. 
(For  example,  when  it  is  done  on  a  TRS-80  Color  Computer,  3  1/2  hours  is  required  to 
generate  an  approximate  boundary  for  a  500-point  image.)  The  required  computation 
time  can  be  reduced  by  including  five  additional  points  at  a  time.  When  the  error  is 
found  to  exceed  UP,  one  point  at  a  time  is  discarded  until  UR  <  UP  .  In  this  way,  the 
computation  time  is  reduced  by  roughly  four.  The  version  of  the  Boundary  Compacter 
listed  in  Appendix  B  is  implemented  thus. 
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CHAPTER  4.  SUMMARY 

This  report  describes  the  algorithms  for  the  detection  and  generation  of  the  binary 
image  boundary  given  the  bit-map  representation  of  the  image  and  for  finding  a  set  of 
polynomial  segments  which  approximate  the  boundary  with  specified  error  tolerance.  In 
the  prototype  automatic  scaling  facility,  these  algorithms  were  implemented  in  BASIC. 
(See  subroutines  Boundary  Eater  and  Boundary  Compacter  in  Appendix  B.) 

4-1.  Performance  of  the  Algorithms 

This  section  describes  two  examples  illustrating  the  performance  of  these 
algorithms.  The  figures  described  here  were  generated  on  a  TRS-80  Color  Computer 
and  an  EPSON  MX-80  Printer.   Appendix  A  describes  the  operating  instructions. 

Figure  4  shows  a  binary  image  which  has  a  variety  of  contours  and  sharp  corners. 
It  was  constructed  from  quadratic  segments.  (The  program  used  to  generate  the  image 
is  also  listed  in  Appendix  B.  Lines  50  -  60  generate  the  tic  marks  on  the  border  at 
intervals  of  20.  Lines  100  -  130  generate  the  various  parabolic  segments.  Line  135 
generates  a  fill  of  the  interior.)  The  bit-map  of  this  image  was  used  as  input  to  the 
Boundary  Eater  Subroutine.  The  Boundary  Eater  subroutine  extracted  the  boundary  of 
the  image  as  shown  in  Figure  7. 

The  subroutine  Boundary  Compacter  then  finds  polynomial  segments  which 
approximate  the  boundary.  Table  1  lists  the  data  generated  as  cubic  segments  are  grown 
to  fit  the  boundary.  In  this  table,  ER  denotes  the  average  absolute  error  on  an  interval 
and  UR  denotes  the  maximum  error  on  an  interval.  The  first  10  lines  list  the  values  of 
ER  and  UR  for  each  group  of  five  boundary  points  considered.  The  maximum  error  was 
found  to  be  smaller  than  the  user  specific  maximum  tolerated  error  (1.3  in  this  case.)  In 
the  11th  line,  the  value  of  UR  exceeded  1.3.  The  following  two  lines  show  the  values  of 
UR  and  ER  as  one  point  at  a  time  was  removed  from  consideration.  The  value  of  UR 
was  again  found  to  be  less  than  1.3  in  line  12.  Line  13  gives  the  values  of  errors  as  a  new 
group  of  five  boundary  points  were  included.  Table  2  summarizes  the  information  on 
the  different  segments  generated  by  the  Boundary  Compacter.  For  example,  the  first 
segment  consists  of  boundary  points  from  0  to  53.  The  coordinates  of  the  starting  and 
ending  points  are  given  by  (X(S),  Y(S))  and  (X(E),  Y(E)),  respectively.  The  coefficients 
of  the  cubic  polynomials  are  given  by  the  values  of  A  and  B  in  the  first  row.  This 
information  is  needed  to  reconstruct  the  entire  boundary.    Specifically,  for  each  segment, 
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Figure  7.  Boundary  Generated  by  Boundary  Eater 
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UR  = 

.379848244 

ER  = 

H  ~  1 C  ■  I  7  C  1 

*  i  .'  i  JO  ■  0  Ji 

ER  = 

.202240307 

UR  = 

.414654261 

ER  = 

.227027074 

UR  = 

.55518995 

ER  = 

.23260026  UR  = 

.552716996 

ER  = 

.240194617 

UR  - 

.594553335 

ER  = 

.208332425 

UR  = 

.566982197 

ER  = 

.272188676 

UR  = 

.943136127 

ER  = 

.241792074 

UR  = 

. 60  761 4 156 

ER  = 

.25970131  UR  = 

.741482634 

ER  = 

.335777036 

UR  - 

.346316725 

ER  = 

.62216693  UR  = 

2.03155914 

ER  = 

.462516072 

UR  = 

1.51663388 

ER  = 

.35062976  UR  = 

1.01514314 

ER  = 

1  Trt"J7--)7C  MP,  _ 

.296437767 

ER  = 

UR  = 

.531935764 

ER  = 

. 22246 j963 

UR  = 

.617656428 

ER  = 

.186466821 

UR  = 

.635330388 

ER  = 

.225501594 

UR  = 

.773245005 

ER  = 

.227226287 

UR  = 

.  /  11  ilJi. 

ER  = 

,226065273 

UR  = 

.812536391 

ER  = 

.233905326 

UR  = 

.770322391 

ER  = 

.234637308 

UR  = 

. 770273026 

ER  = 

.248309924 

UR  = 

.727645327 

ER  = 

.261717261 

UR  = 

.340892526 

ER  = 

.261732355 

UR  = 

.709957583 

ER  = 

.247729389 

UR  = 

.715134679 

ER  = 

.263615909 

UR  = 

.708693612 

ER  = 

.271316683 

UR  = 

.72633242 

ER  = 

.312760945 

UR  = 

1 .2155355 

ER  = 

.314670033 

UR  = 

.997031434 

ER  = 

.864626712 

UR  = 

3.97035329 

ER  = 

.639550961 

UR  = 

2.79073362 

ER  = 

.445393187 

UR  = 

1.70846895 

ER  = 

.316542415 

UR  = 

.932250155 

ER  = 

0                UR  = 

0 

ER  = 

.142788436 

UR  = 

.449412781 

ER  = 

.17675141  UR  = 

.510605028 

ER  = 

.191724772 

UR  = 

.539347639 

ER  = 

.194349524 

UR  = 

.550827266 

ER  = 

.227907526 

UR  = 

.750169185 

ER  = 

.292136233 

UR  = 

1.38110018 

ER  = 

.267401273 

UR  = 

.802060224 

ER  = 

0                UR  = 

0 

ER  = 

.230395767 

UR  = 

.453745874 

ER  = 

.217920354 

UR  = 

.75221239 

ER  = 

.235263434 

UR  = 

.713433005 

ER  = 

.213990781 

UR  = 

.738730705 

Table  1  (Part  1) 
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ER  = 

.23461789  UR  = 

.701202642 

ER  = 

.288880457 

UR  = 

1.09636243 

ER  = 

.279183134 

UR  = 

.869242959 

ER  = 

.306568527 

UR  = 

.954480032 

ER  = 

.282733171 

UR  = 

.842906589 

ER  = 

.345000952 

UR  = 

.978749789 

ER  = 

.38280489  UR  = 

1.2750774 

ER  = 

.381860283 

UR  = 

1.44895442 

ER  = 

.376986516 

UR  = 

1.40506563 

ER  = 

.372817638 

UR  = 

1.3661353 

ER  = 

.370719588 

UR  = 

1.33189632 

ER  = 

.372748895 

UR  = 

1.30183026 

ER  = 

.38280489  UR  = 

1.2750774 

ER  = 

0      UR  = 

0 

ER  = 

.163172984 

UR  = 

.489238558 

ER  = 

.216970337 

UR  = 

.668754398 

ER  = 

.24990067  UR  = 

.706899183 

ER  = 

.204722329 

UR  = 

.524873336 

ER  = 

.203931045 

UR  = 

.514338793 

ER  = 

.236984343 

UR  = 

.648389815 

ER  = 

.237156489 

UR  = 

.663141601 

ER  = 

.244033019 

UR  = 

.759129206 

ER  - 

.248349713 

UR  = 

.816529769 

ER  = 

.246420879 

UR  = 

.850575885 

ER  = 

.295914282 

UR  = 

.981974668 

ER  = 

.271739502 

UR  = 

1.0327295 

ER  = 

.284193043 

UR  = 

1.09712407 

ER  = 

.324817926 

UR  = 

1.24394043 

ER  = 

.425624298 

UR  = 

1.43813997 

ER  = 

.389152794 

UR  = 

1.37799632 

ER  ' 

.448938607 

UR  = 

1.41807838 

ER  = 

.405981272 

UR  = 

1.3593946 

ER  = 

.361617136 

UR  = 

1.30098303 

ER  = 

.324817926 

UR  = 

1.24394043 

ER  = 

.110849409 

UR  = 

.284433906 

ER  = 

.118940079 

UR  = 

.309980591 

ER  = 

.205964267 

UR  ■ 

.478622288 

ER  = 

.221155891 

UR  = 

.582018434 

ER  = 

.190380066 

UR  ■ 

.493356459 

ER  = 

.164895437 

UR  ■ 

.489725185 

ER  = 

.16138711  UR  - 

.494159223 

ER  = 

.172170238 

UR  * 

.487455033 

ER  = 

.17095724  UR  » 

.481125895 

ER  ■ 

.193442374 

UR  * 

.534429748 

ER  * 

.196603322 

UR  « 

.519499039 

ER  » 

.290660989 

UR  ■ 

1.06120439 

ER  * 

.347156175 

UR  * 

1.17646111 

ER  ■ 

.641975514 

UR  - 

1.91005969 

Table  1  (Part  2) 
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ER  = 

.648248797 

UR  = 

1.95247873 

ER  = 

.416794277 

UR  = 

1.30399358 

ER  = 

.399382956 

UR  = 

1.28412306 

ER  = 

.177567352 

UR  = 

.379848244 

ER  = 

.168379723 

UR  = 

.430910956 

ER  = 

.177552679 

UR  = 

.485034066 

ER  = 

.205525705 

UR  = 

.443171589 

ER  = 

.212508105 

UR  = 

.44465278 

ER  = 

.236985102 

UR  = 

.707095407 

ER  = 

.251020596 

UR  = 

.687387859 

ER  = 

.273676504 

UR  = 

.801894931 

ER  = 

. 306206965 

UR  = 

.932947065 

ER  = 

.33593672  UR  = 

1.04285935 

ER  = 

.404970467 

UR  = 

1.30154677 

ER  = 

.33741846  UR  = 

1.2513366 

ER  = 

0      UR  = 

0 

ER  = 

.142788436 

UR  = 

.449412781 

ER  = 

.137865537 

UR  = 

.411234053 

ER  = 

.157016348 

UR  = 

.48707989 

ER  = 

.16960535  UR  = 

.464473602 

ER  = 

.156268968 

UR  = 

.521708434 

ER  = 

.135319249 

UR  = 

.555308994 

ER  = 

.193304898 

UR  = 

.566324501 

ER  = 

.20342355  UR  = 

.534303207 

ER  = 

.19958201  UR  = 

.530552704 

ER  = 

.210415783 

UR  = 

.531869313 

ER  = 

.247255424 

UR  = 

.784290909 

ER  = 

.22153644  UR  = 

.613647949 

ER  = 

.21827324  UR  = 

.636992913 

ER  = 

.213903011 

UR  = 

.631753831 

ER  = 

.225543026 

UR  = 

.392507443 

ER  = 

.325420955 

UR  = 

1.49247134 

ER  = 

.266040915 

UR  = 

1.07863389 

ER  = 

.177567352 

UR  = 

.379848244 

ER  = 

.208770874 

UR  = 

.393110046 

ER  = 

.172032454 

UR  = 

.436345162 

ER  = 

.22942356  UR  = 

.430233006 

ER  = 

.191178846 

UR  = 

.536884233 

ER  = 

.727073653 

UR  = 

2.99100406 

ER  = 

.541357382 

UR  = 

2.31915107 

ER  = 

.363116314 

UR  = 

1.56i!67238 

ER  = 

.239194894 

UR  = 

1.04098208 

TIMER  VALUE  =  4587 

Table  1  (Part  3) 
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only   the   coordinates  of  the  starting  point,   X(S),  Y(S),  and   the  parameters  of  the 
polynomial  are  retained. 

The  approximate  boundary  reconstructed  from  the  data  in  Table  2  is  shown  in 
Figure  8.  This  approximate  boundary  is  generated  from  9  segments  and  a  maximum 
tolerated  error  of  1.3.  Figure  9  shows  a  different  approximate  boundary  with  a  smaller 
maximum  tolerated  error  (0.7).  The  17  segments  of  this  approximate  boundary  are  listed 
in  Table  3.    In  this  example,  both  approximate  boundaries  are  visually  acceptable. 

Another  example  considered  was  the  binary  image  whose  boundary  is  shown  in 
Figure  10a.  This  boundary  consists  of  ellipic  and  circular  segments  and  was  generated 
from  the  data  in  line  5  of  the  program  listed  in  Appendix  B.  The  resultant  boundary 
found  by  Boundary  Eater  is  shown  in  Figure  10b.  With  maximum  tolerated  error  of  0.9, 
the  approximate  boundary  is  given  by  the  data  in  Figure.  11.  In  this  case,  we  can  see  an 
undesirable  result  where  the  last  segment  wraps  around  and  covers  part  of  the  first.  The 
resultant  approximate  boundary  is  not  satisfactory  visually  because  the  symmetry  in  the 
original  image  was  not  preserved.  Figures  12-15  show  the  approximate  boundaries  of 
this  image  when  the  maximum  tolerated  error  was  0.5,  0.3,  0.15  and  0.0001,  respectively 
Unfortunately,  the  quality  of  the  approximate  boundary  suffers  severely  when  the 
maximum  tolerated  error  is  too  small.  This  may  be  due  to  problems  in  the  subroutine 
Better  Plotter  which  graphed  the  segments. 

4-2.   Further  Work 

To  scale  a  binary  image,  the  approximate  image  boundary  found  by  the  Boundary 
Compactor  needs  to  be  reduced  to  the  desired  size.  The  smaller  image  at  the  desired 
resolution  may  then  be  obtained  by  inking  the  interior  of  the  reduced  boundary.  The 
codes  required  to  carry  these  two  steps  are  not  included  in  this  report.  Several 
experiments  were  carried  to  determine  the  visual  qualities  of  reduced  images  produced  in 
this  fashion.  When  the  numbers  of  pixels  in  all  portions  of  an  image  are  large,  the 
quality  of  the  reduced  image  in  coarser  resolutions  is  usually  good  as  expected.  When 
some  parts  of  the  image  contain  very  few  pixels  (e.  g.,  around  5), the  quality  of  the 
resultant  image  was  found  poor  in  some  cases.  The  poor  quality  is  due  to  "quantization 
noise"  in  these  cases.  In  order  to  assure  high  visual  quality  in  the  final  result,  supports 
which  allow  the  user  to  edit  the  bit-map  representation  generated  automatically  by  the 
scaling  facility  must  be  provided.  These  supports  need  to  be  implemented  in  the  future. 
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Figure  8.   Approximate  Boundary  with  Maximun  Tolerated  Error  =  1.3. 


24 


I 


\ 


nH 

f 

i 

■ 

_■ 

I 

,-' 

i 

\ 

/ 

i 
i 

«, 

f 

1 

\ 

/ 

■ 

H 

\. 

a 
ii 

1 

n 

■-« 

u 
tat 
■i 

) 

\ 


\ 

u 


•■Si 

nan  » 

aa  aan 


I 


/ 


/ 

i 

i 

/ 


Figure  9.  Approximate  Boundary  with  Maximum  Tolerated  Error  =  0.7 
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Figure  10a.   Actual  Boundary  of  a  Symmetrical  Image 
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Figure  10b.   Boundary  of  the  Image  in  10a  Found  by  Boundary  Eater 
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Figure  13.   Approximate  Boundary  with  Maximum  Tolerated  Error  =  0.3 
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Figure  14.   Approximate  Boundary  with  Maximum  Tolerated  Error  =  0.15 
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APPENDIX  A.  OPERATIONS   OF  ASF 

The  programs  listed  in  Appendix  B  were  ran  on  a  TRS-80  Color  Computer.  An 
EPSON  MX-80  Printer  was  used  to  plot  figures  described  in  Chapter  4.  The  operating 
instructions  are  listed  below.  Lines  of  text  in  italics  indicate  commands  or  responses 
entered  by  the  user.  Similarly,  actions  taken  by  the  user  are  shown  in  italics.  Bold  face 
text  indicates  prompts  from  the  system.  The  text  in  Roman  indicates  comments. 

(1)  CLEAR   100,  32500 

This  command  is  entered  before  loading  the  program,    clearing  space  for  a 
machine  language  subroutine  —  see  (7). 

(2)  PMODE  4,1:  PCLS  5:  RUN 

This  command  sets   proper  graphics   mode,  clears  graphics  screen,  runs  the 
program. 

(3)  MAGNIFICATION?  1. 

(4)  Circle  is  drawn,  BREAK 

This  action  returns  the  program  to  user  control. 

(5)  PCLS  5:   GOSUB  1100 

(6)  The  image  shown  in  Figure  4  *»  drawn,  BREAK 

(7)  GOSUB  18000 

This  command  loads  a  machine  language  routine  which  assists  subroutine  12000 
in  posting  the  screen  to  the  printer. 

(8)  GOSUB  12000 

This  command  posts  the  screen  to  the  printer. 

(9)  GOSUB  10000 

This  program  finds  the  boundary  of  the  screen  image.  The  user  is  prompted  to 
enter  data  needed.  The  data  entered  here  include  the  coordinates  of  the  starting 
point  (a  point  on  the  boundary),  the  direction  of  the  associated  boundary  vector, 
and  the  relationship  between  the  boundary  point  and  the  boundary  vector. 
START  PT?  83,40 

83  and  40  are  coordinates  of  boundary  point. 
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START  VECTOR  DIR   U,  D,  L,  R,  ?  R 

U,  D,  L,  and  R  means  Up,  Down,  Left,  or  Right,  respectively. 
START  PT-VEC  REL  G,  S?  G 

G  and  S  means  Greater  or  Smaller  than  the  vector,  respectively. 

This  subroutine  when  completed,  erases  the  original  image  and  draws  the 
boundary  on  the  screen.  The  Coordinates  of  the  boundary  points  are  saved  in 
arrays  X  and  Y,  each  with  TP  elements. 

(10)  GOSUB  12000 

This  command  causes  the  image  boundary  to  be  dumped  to  the  printer. 

(11)  GOSUB  21000 

This    command    starts    the    Boundary    compacter    to    find    the    approximate 
boundary.  The  user  is  prompted  to  enter  the  maximum  tolerated  error. 
STARTING  PT  #?  0 
maximum  tolerated  error?  1.3 

A  value  between  .25  and  1.9  seems  reasonable.  As  this  program  runs,  error 
information  on  various  segments  is  printed.  When  this  subroutine  is  completed,  the 
information  about  the  curves  which  approximate  the  boundary  is  erased  and  the 
approximate  boundary  is  graphed. 

(12)  GOSUB  12000 

This  command  sends  the  approximate  boundary  to  the  printer. 

(13)  By  comparing  the  output  generated  in  step  No.  12,  the  approximate  boundary,  to 
the  output  generated  in  step  No.  10,  the  original  boundary,  the  user  can  decide  if 
he  is  satisfied  visually  by  the  approximation.  If  he  is  not  satisfied,  he  may  repeat 
step  No.  11  with  a  smaller  tolerated  error. 
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APPENDIX  B.  PROGRAMS  USED  IN  ASF 


1  REM  COPYRIGHT  ©HAY  7,1983  BY  JEFFREY  S.  ELY 

2  INPUT"MA6NIFICATI0NH;MG 
5  DATA  C, 100, 100,20,. 

20  DATA  1,50,20,110,20 
22  DATA  L, -,-,70,60 
24  DATA  L,-, -,130,60 
26  DATA  L,-,-, 170, 20 
28  DATA  L,-,-, 230, 20 
30  DATA  L,-,-, 130,70 

32  DATA  L, -,-,20,70 
34  DATA  L,-,-,20,50 
36  DATA  L,-, -,50,20 

33  DATA  F, 60,40,. 

50  DATA  P,0,0,P,20,0,P,40,0,P,60,0,P,80,0,P,100,0,P,120,0,P,140,0,P,160,0,P,180,0,P,200,0,P,220,0,P,240,0 

60  DATA  P,0,20,P,0,40,P,0,60,P,0,80,P,0,100,P,0,120,P,0,140,P,0,160,P,0,180 

100  DATA  0,83,40,-1,0,27,95,0,1 

105  DATA  Q, -,-,-,-,117, 137, 9, -12 

110  DATA  9, -,-,-,-, 176, 98,4,1 

115  DATA  Q,-, -,-,-,228,70,-2, -15 

120  DATA  0,228,70,-1,0,210,43,0,1 

125  DATA  0,210,43,-70,-43,140,60,-55,50 

130  DATA  0,-,-,-,-, 83, 40, -6,-19 

135  DATA  F,60,80,. 

930  mmmmmmmmmmmmmmmmmmm$ 

950  REH  MAIN  PROGRAM  TO  READ  DATA  AND  INTERPRET  COMMANDS 

970  mimwmfflmmmminnnmmmMmimmm 

1000  SN=5:DIM  SY$(SN) 

1020  SY$(1)="L":SY$(2)="C':SY$(3)="P":SY*(4)="FH:SY$(5)="S" 

1030  LM=550:DIM  X(LM),V<LH):REH  USED  BY  BOUNDARY  EATER  TO  STORE  POINTS  EATEN 

1045  LN=45:DIH  S(LN),E<LN),A(LN),B(LN) 

1100  PMODE  4,1: SCREEN  1,1 

1200  READ  BS* 

1220  IF  DS$="."  THEN  1400:REM  END  OF  COMMAND 

1240  SI=0 

1260  FOR  K=l  TO  SN 

1280  :   IF  SY$«K)=DS*  THEN  SI=K 

1300  NEXT  K 

1320  IF  SI-0  THEN  PRINT  "ILLEGAL  SYMBOL ";DS$: GOTO  1500 

1340  ON  SI  GOSUB  4000,5000,6000,7000,8000 

1360  GOTO  1200 

1400  REM   13000:60SUB  12000:REM  POST  SCREEN  TO  PRINTER 

1450  GOTO  1450 

1500  STOP 


3930  REMMM3MMMMi»MMMMMWMMMMMMMMMJMittM 

3950  REN  LINE 

3970  i£.mimfflmmmm%mtfflmmmm3mmmm 

4000  READ  LW*,LX*,LY,L2 
4002  LY=LYtH6:LZ=LZMJG 

4005  IF  LW$="-"  THEN  L«=QX  ELSE  LH=VAL!Ltt*>tHG 
4010  IF  LS*="-"  THEN  LX=9Y  ELSE  LX=VAULXt>tH6 
4020  LINE<LW,LX}-!LY,LZ>, PRESET 
4040  QX=LY:BY=LZ 

4060  DX=LY-LW:DY=LZ-LX 

4100  RETURN 

4930  RErtiii335333WS35wJJ33iJSSS3^i5i35#ij>3i**ii3:iW^3ii*5id 

4950  REM  CIRCLE 

4970  REt1^iiii333W3£*3*J353i'5Wi*iviW-i5i^i3w^JMi»i5ii35i.}J 

5000  READ  CA,CB,CR 

5015  CA=CAiM6:CB=CB*M6:CR=CRtMe 

5020  CIRCLE <CA,CBi,CR,0 

5100  RETURN 

5930  REH389393ft933ttMtt3tf$ft3333d9d3?3MSdd33d3M333i3d9333 

5950  REM  POINT 

5970  mmmmmmmmmmmmmmmmmmm 

6000  READ  PX,PY 

6015  PX=PXtM6:PY=PYtM6 

6020  PSET(PX,PY,0) 

6100  RETURN 

6930  Rm-mmmmmm$mmmmmfaimmwdjMim 

6950  REM  FILL 

6970  REMii3WJJ53J3*3J33iS;J33i.il5wvi535ii33SiliW5J5wi5Wi55i* 

7000  READ  FX,FY 

7015  FX=FX*M6:FY=FY*H6 

7020  PAINT (FX, FY), 0,0 

7100  RETURN 

7930  REHMWiWMiSiMWMSSdMwiJSMWWSMWSWMSiSMSWM 

7950  REM  QUADRATIC 

7970  mmmmmmmmmmmmzmmmmmMm 

8000  READ  Xl*,Yl*,Ui$,Vi*fX2,Y2,U2,V2 

8001  X2=X2tHG:Y2=Y2tMG 

8002  IF  Xl*="-"  THEN  X1=QX  ELSE  Xl=VftL(Xl$)W8 
8004  IF  Ylt="-"  THEN  Y1=QY  ELSE  Y1=VAL(Y1$.UMB 
8006  IF  Ult="-"  THEN  U1=DX  ELSE  Ul=VAL(Ult) 
8008  IF  Vll="-"  THEN  V1=DY  ELSE  Vl=VAL(Vlt> 
8010  REM  DETERMINE  COEFFICIENTS  A,B,C,D 

8020  B=X2-X1;D=Y2-Y1 

8040  Ti=VUB-UltD  :  T2=V2tB-U2»D 

8060,DI=VltU2-V2tUl 

8080  IF  DI=0  THEN  PRINT  "DI=0n:80T0  1500 

8300  C=(V2»T1+VUT2)/DI 

8120  A=(U2tTl+UltT2)/DI 

8130  REM  DETERMINE  LARGEST  PARAMETER  SPACING, H, FOR  CONTIN  PLOT 

8140  G1=ABS(B-A):G2=ABS(B+A):G3=ABS!D-C):G4=ABS(D+C) 

8160  HX=G1 

8180  IF  62>MX  THEN  MX=G2 
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8200  IF  63>HX  THEN  HX=S3 

8220  IF  64>HX  THEN  MX=S4 

8240  IF  MX=0  THEN  PRINT  "HX=0":S0T0  1500 

8260  H=1.9/MX 

8270  REH  PLOT  THE  CURVE 

8273  X4=INT(X2+.5):Y4=INT(Y2+.5>:PSET(X4,Y4,0} 

8275  X3=INT(XH-.5):Y3=INT(Y1+.5):LINE(X3,Y3)-(X3,Y3),PRESET 

8280  T=H 

8285  IF  T>=1  THEN  8335 

8300  :   Rl=T*iT-l> 

8320  :   Pl=A*Rl+BtT+Xl:PI=INTiPH-.5}:XD=flt(T+T-l)+B 

8340  :      P2=C»Rl+DtT+Yl:PJ=I»ITtP2+.5hYD=Ct<W-l>+D 

8343  IF  (INT(P1)=L1)  AND  (INTIP2>=L2>  THEN  8380 

8345  :      L1=INT(P1):L2=INT<P2}:H1=L1+1:H2=L2<-1 

8350  :   IF  ABSiVD) >ABSiXD)  THEN  3363 

8353  :   9U=INT(P2+(YD/XD>*Ul-Pl)+.5) 

8354  LINE-iLl,8U>, PRESET 
3360  :   GOTO  8380 

8363  REH  60  ON 

3366  :   B»=INT(Pl+<XB/YD}»<L2-P2>+.5) 

8367  LINE-(BW,L2>, PRESET 

8330  T=T+H:60T0  8285 

8385  LINE-(X4,Y4), PRESET 

8390  9X=X2:BY=Y2 

3395  BX=U2:DY=V2 

8400  RETURN 

9930  mmmmmmmmdmmdiJim^*mi*ummM 

9950  REM  BOUNDARY  EATER 

9970  mmmmmmmmmmmmm^immmmM 

10000  INPUT  "START  PT";X,Y 

10005  INPUT  "START  VECTOR  DIR  U,D,L,R";DRt 

10010  INPUT  "START  PT-VEC  REL  G,S";RL$ 

10020  T6=.5:T7=.5:T8=.5:T9=.5 

10022  IF  RL$="S"  THEN  T6=-.5:T7=-.5:T8=-.5:T9=-.5 

10024  IF  DR*="R"  THEN  T6=-.5:T3=.5 

10026  IF  DRt="D"  THEN  T7=-.5:T9=.5 

10028  IF  DRt="L"  THEN  T8=-.5:T6=.5 

10030  IF  DRt="U"  THEN  T9=-.5:T7=.5 

10035  BX=X+T6:8Y=Y+T7;EX=X+T8:EY=Y+T9 

10050  WH=5:BL=0:REH  WHITE, BLACK 

10060  TP=-1:REH  INIT  POINTER  TO  TOP  OF  POINT  STORE 

10070  TP=TP+1:X<TP)=X:Y(TP)=Y 

10090  REH  BE6IN  LOOP-FIND  NEXT  PDINTtttttttttUt 

10100  REH  El  DEPENDS  ONLY  ON  DR$,INIT  TO  0 

10110  E1=0 

10120  IF  DR*="L"  THEN  El=-1 

10130  IF  DR$="R"  THEN  El=l 

10140  REH  E2  DEPENDS  ONLY  ON  DR$,INIT  TO  0 

10150  E2=0 

10160  IF  DR*="U"  THEN  E2=-l 

10170  IF  DR$="D"  THEN  E2=l 
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10190  Dl=l 

10200  IF  DR*="L"  THEN  Dl=-1 

10210  IF  RL«="6"  AND  DR$<>"R"  THEN  Dl=-1 

10230  D2=l 

10240  IF  DR$="U"  THEN  D2=-l 

10250  IF  RL*="S"  AND  DRJO°D"  THEN  D2=-l 

10270  REM  COMPUTE  DO  AND  OT 

10280  IF  PP0INT<X+D1,Y+D2)=BL  THEN  D0=1  ELSE  DQ=0 

10290  IF  PP0INT(X+E1,Y+E2)=SL  THEN  0T=1  ELSE  0T=0 

10310  REM  COMPUTE  W1,W2 

10320  W1=0:W2=0 

10330  IF  D0=1  THEN  W1=D1:W2=D2 

10340  IF  D0=0  AND  0T=3  THEN  W1=E1:W2=E2 

10350  REM  COMPUTE  VI, V2 

10360  IF  RL$="S"  THEN  J6=-1:J7=1:J8=-1:J9=1 

10370  IF  RLJ="6"  THEN  J6=l:J7=-l:J3=l:J9=-i 

10330  IF  DR$="L"  OR  DR$="R"  THEN  J6=0:J7=0 

10390  IF  DR$="U"  OR  DR*="D"  THEN  JS=0:J9=0 

10420  V1=E1:V2=E2 

10430  IF  D0=1  THEN  V1=J7:V2=J9 

10440  IF  D0=0  AND  0T=0  THEN  V1=J6:V2=J8 

10450  REM  COMPUTE  NEW  DR$ 

10455  NR$=DR$:REM  SAVE  OLD  SO  NEW  RL$  CAN  USE  IT  BELOW 

10460  IF  Vl=-1  THEN  DR$='*L" 

10470  IF  Vl=l  THEN  DR$="R" 

10480  IF  V 1 =0  AND  V2=-l  THEN  Dfit="Ll" 

10490  IF  V1=0  AND  V2=l  THEN  DR$="DU 

10500  REM  COMPUTE  NEW  RL$ 

10510  ja$=u3u:J9$="B" 

10520  IF  NR$="U"  OR  NRt="L*  THEN  J8*="G":J9t=aS" 

10530  IF  D0=1  THEN  RL$=J9$ 

10540  IF  D0=0  AND  0T=0  THEN  RL$=J8« 

10600  REM  UPDATE  POINT  i  VECTOR  END 

10610  X=X+«i:Y=Y+W2:E)(=EX+Vi:E¥=EYn,2 

10620  REM  STORE  POINT  IF  NEW 

10630  IF  IXO)((TP))  OR  (YOYITP))  THEN  TP=TP+1:X(TP)=)!:Y('TP)=Y 

10640  rem  test  for  ENDmttmmmtmmtm 

10650  IF  (EXOBX)  OR  (EYOBY)  THEN  10110 

10660  PCLS  5 

10664  FOR  K=0  TO  240  STEP  20:PSET(K,0,BL):NEn  K 

10666  FOR  K=0  TO  180  STEP  20:PSET<0,K,BL>:NE)iT  K 

10668  SCREEN  1,1 

10670  FOR  K=0  TO  TP:PSET (X (K) , Y<K> ,BL) :NEXT  K 

10700  RETURN 

1 1 930  REMMM3MMS3 JS mi immmmm mmmm mi mi 

11950  REM  POST  SCREEN  TO  PRINTER 

11970  re  mm  mmm i  a mmummm m am m mmi  n 

12000  PRINT#-2,CHRt(27);"A";CHR$(0);CHR$«15) 
12020  AD=1536 
12040  FOR  R=l  TO  64 
12060  FOR  C=i  TO  32 
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12070  Bl=PEEK(AD>:B2=PEEK<AB+32):B3=PEEKlAD+64> 

12080  POKE  32500, Bl: POKE  32501, B2:P0KE  32502, B3 

12090  IFiBl  AND  B2  AND  B3)=255  THEN  PRINT#-2,"    °;:GOTO  12180 

12100  FOR  U=l  TO  4:EXEC  32504:PRINTi-2.CHRt(PEEK(32503)*160);:NEXT  U 

12180  AD=AD+1 

12200  NEXT  C 

12220  PRINT#-2,"  " 

12230  AD=AD+64 

12240  NEXT  R 

12250  RETURN 

12930  mmmmmmmmmmmmmmmmmma 

12950  REM  HACH  LAN6  SCREEN  DUMP  SUPPORT 

1 2970  REMiiiiiiSii h Ji m s mn namm mmm i w nmm 

13000  DATA  121,126,244,118,126,247 

13010  DATA  121,126,244,118,126,247 

13020  DATA  121,126,245,118,126,247 

13030  DATA  121,126,245,118,126,247 

13040  DATA  121,126,246,118,126,247 

13050  DATA  121,126,246,118,126,247 

13060  DATA  67,118,126,247 

13070  DATA  67,118,126,247 

13080  DATA  115,126,247 

13090  DATA  57 

13100  FOR  1=32504  TO  32551 

13110  READ  A 

13120  POKE  I, A 

13130  NEXT  I 

13200  RETURN 

1 7930  mmmmmnmnmimmmmmmmmmt 

17950  REM  LEAST  SQUARES  CUBIC  C0EFFS 

1 7970  mmmmmmmmmmmimmmiimmij 

18000  INPUT  "S,E";S,E 

18010  Bl=X(E)-X<S):B2=Y<E)-Y(S):BB=S0R(BltBl+B2»B2> 

1 8020  S1=0: 32=0: S3=0: S4=0: S5=0 

18030  K=S 

18040  REH  BEGIN  LOOP  (D0NT  DELETE) 

18050  A1=X(K)-X(S):A2=Y<K)-Y<S) 

13060  XB=IAl*Bl+A2tB2)/BB:YB=(Al»B2-Bl»A2)/BB 

18070  XD=XB-BB:Q0=XB»XD 

18080  Ql=60tYB:Sl=Sl+91 

18090  Q2=Q1*XB:S2=S2+Q2 

18100  Q3=Q0*Q0:S3=S3+fl3 

18110  Q4=Q3tXB:S4=S4+Q4 

18120  Q5=fi4tXB:S5=S5+B5 

13130  IF  K=E  THEN  18155 

18140  K=K+1:IF  K>TP  THEN  K=0 

18150  GOTO  18040 

18155  REH  GO  ON  (D0NT  DELETE) 

18160  DT=S3tS5-S4tS4 

18170  IF  DT=0  THEN  PRINT  'NO  SOL  TO  LEAST  SQUARES":60TQ  18200 

18180  A=(32tS3-SUS4)/DT:B=(SltS5-S2tS4)/DT 

18200  RETURN 
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18330  9£mmm*mwwnwwwnmmmm9mMnm 

18350  REH  AVE  ABSOLUTE  AND  HA)!  ERROR  FOR  LEAST  SQ  CUBIC 

1 83?o  re ii id m w a m i j m n a m i » m S3 nm m mm m 

18400  ER=0:UR=0:CN=0 

18410  K=S 

18420  REH  LOOP  (DONT  DELETE) 

18430  CN=CN+1 

18440  A1=K(K}-X(S)jA2=Y(K>-Y<S5 

18450  XB=(AltBl*A2tB2)/BB:YB=<AUB2-Bl*A2)/B8 

18460  YT=  ( AtXB+B)  t)£B*  t XB-BB) 

13470  LR=ABS(YT-YB> 

13480  ER=ER+LR 

13490  IF  LR>UR  THEN  UR=LR 

13500  IF  K=E  THEN  18530 

18510  K=K+1:IF  IOTP  THEN  K=0 

18520  GOTO  18420 

13530  REH  60  ON  (DONT  DELETE) 

18540  ER=ER/CN 

18550  PRINTf-2,"ER  =  ";ER,"UR  =  ";UR 

18600  RETURN 

19030  m%M%M%%%%%m%%fflimwmmmnnwmmnwti 

19050  REH  BETTER  PLOTTER  FOR  LEAST  SQUARES  CUBIC 

19070  REHi3SiWWSJSWJi35WWW53W3WW^i3iJMWJ3iwSi55i 

19100  IF  A=0  AND  B=0  THEN  LINE(SUS),Y(S))-()i(E},Y(E)i,PRESET:GOTO  19450 

19105  IF  A=0  THEN  81=0:60T0  19140 

19110  X9=(A*BB-B)/(3tA) 

19120  IF  (X9>BB>  OR  (X9<0>  THEN  81=0:6010  19140 

19130  61=ABS{3tAtX9tX9+2t(B-AtBB)«X9-BtB8) 

19140  REH  60  ON  (DONT  DELETE) 

19150  62=ABS(BtBB) 

19160  63=ABS((At3B+B)tBB) 

19170  HX=61 

19180  IF  62>HX  THEN  HX=62 

19190  IF  B3>HX  THEN  WX=63 

19200  IF  MX=0  THEN  PRINT  "HX=0":60T0  1500 

19210  H=1.9/HX:REH  MAX  XBAR  CHAN6E  FOR  600D  GRAPH 

19220  REH  PLOT  THE  CURVE 

19230  X4=INT(X(E)+.5):Y4=INW(E}+.5)jPSETU4,Y4,0) 

19240  X3=INT(X(S>  +  . 5) :Y3=INT<Y<S>«-.5>:LINEIX3,Y3>-IX3,Y3>, PRESET 

19250  XT=H:Z1=X3jZ2=Y3:L1=X3jL2=Y3 

19260  IF  XT>=BB  THEN  19410 

19270  YT=(AW+B)Wt()(T-BB) 

19230  Pl=X(S)t(YTtB2+XT»Bl)/BB:PI=INT(Pl+.5):XD=Pl-Zl;Zl=Pl 

19290  P2=Y<S)+(XT»B2-YTtBl)/BB:PJ=INT(P2+.5):YD=P2-22:Z2=P2 

19300  IF  (INT(P1)=L1)  AND  (INT(P2)=L2)  THEN  19390 

19310  L1=INT(P1>:L2=INT(P2):H1=L1+1:H2=L2+1 

19320  IF  ABS!YD))ABS(XD)  THEN  19360 

19330  QU=INT(P2+(YD/XD>*(L1-P1)0) 

19340  LINE-(L1,BU), PRESET 

19350  GOTO  19390 

19360  REM  60  ON  (DONT  DELETE) 

19370  BW=INT(PH-(XD/YD)t(L2-P2>+.5) 
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19380  ljne-(b»,l2>, preset 
19390  xt=xt+h 

19400  SOTO  19260 

19410  LINE-(X4,Y4), PRESET 

19450  RETURN 

20930  REMM3Mtt«M«MittMWM3iMMa«MM3MMM«aMM 

20950  REM  FASTER  BOUNDARY  COMPACTOR 

20970  REHtt9W339WW33attdWWftftft33d393ttWW33WWW3 

21000  INPUT  "STARTING  PT  i";SR 

21003  INPUT  "MAX  ERROR  TOLERATED"; UP 

21006  JS=5:REM  JUMP  SIZE 

21010  IP=-1:CR=0:S=SR 

21015  TIMER=0 

21020  REM  REPEAT  FINDING  THE  LARGEST  SEGMENTS  (DONT  DELETE) 

21040  E=S 

21050  REM  REPEAT  JUMPING  UNTIL  TOO  BIG  (DONT  DELETE) 

21060  CR=CR+JS:E=E+JS:IF  E>TP  THEN  E=E-TP-1 

21070  GOSUB  18010:G0SUB  18400:REM  OBTAIN  A,B,UR,ER 

21030  IF  UR<UP  THEN  21050 

21090  REM  UR  IS  NOW  TOO  BIG  , START  STEPPING  DOWN  (DONT  DELETE) 

21100  CR=CR-1:E=E-1:IF  E<0  THEN  E=TP 

21110  GOSUB  18010:G0SUB  18400 

21120  IF  UR>UP  THEN  21090 

21130  IP=IP+1:S(IP)=SjE(IP)=E:A(IP)=A:B(IP)=B:REM  SAVE  SEGMENT  INFO 

21140  S=E 

21150  IF  CR<TP  THEN  21020 

21155  PRINTt-2, "TIMER  VALUE  =  "; TIMER 

21160  PRINTI-2,"*  OF  POINTS  IN  FIGURE  =  ";TP+1 

21163  PRINT#-2,"NA)(  ERROR  =  ";UP 

21165  PRINTi-2, "SEGMENT  INFO:" 

21170  PRINTi-2, "S',"E","X(S)","Y(S)",")((E)","Y(E)", "A",  "B":PRINTt-2,"  " 

21180  FOR  K=0  TO  IP:PRINT#-2,S(K),E(K),)({S(K)),Y!S(K)),X(E!K)),Y(E(K)),A(K),B(K):NE.U  K 

21200  REM  GRAPH  THE  RESULTS 

21210  SCREEN  1,1:PCLS5 

21220  FOR  K=0  TO  IP 

21230  S=S(K):E=E(K):A=A<K):B=B(K) 

21240  Bl=X(E)-XIS):B2=Y(E)-Y<S):BB=SQR<BltBl+82*B2) 

21250  GOSUB  19100 

21260  NEXT  K 

21300  RETURN 
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